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1. Introduction 

A lattice polygon is a model of ring polymer in a good solvent, and is a useful in the 
examination of the entropy properties of ring polymers in dilute solution [T^ [5] . Ring 
polymers can be knotted [131 121 HO] , and this topological property can be modeled by 
examining knotted polygons in a three dimensional lattice. The effects of knotting (and 
linking) on the entropic properties of knotted ring polymers remains little understood, 
apart from empirical data collected via experimentation on knotted ring polymers 
(for example, knotted DNA molecules [17]) or by numerical simulations of models 
of knotted ring polymers [HI [3]). Knots in polymers are generally thought to have 
effects on both the physical [IH] and thermodynamic properties [35] of the polymer, 
but these effects are difficult to understand in part because the different knot types in 
the polymer may have different properties. 

A polygon in a regular lattice L is composed of a sequence of distinct vertices 
{uq, ai, . . . , a„} such that ajUj+i and a„ao are lattice edges for each j = 0, 1, . . . , n— 1. 
Two polygons are said to be equivalent if the first is translationally equivalent to 
the second. Such equivalence classes of polygons are unrooted, and we abuse this 
terminology by referring to these equivalence classes as (lattice) polygons. In figure fl] 
we display three polygons in regular cubic lattices. The polygon on the left is a lattice 
trefoil knot in the simple cubic (SC) lattice. In the middle a lattice trefoil is displayed 
in the face-centered cubic (FCC) lattice, while in the right hand panel an example of 
a lattice trefoil in the body-centered cubic (BCC) lattice is illustrated. 
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Figure 1. Lattice trefoils (of knot type 3^ in the standard knot tables 6 ) in 
three dimensional cubic lattices. On the left a lattice trefoil is embedded in the 
simple cubic lattice. In the middle a lattice trefoil in the face-centered cubic 
lattice is illustrated, while the right panel is a realisation of a lattice trefoil in the 
body-centered cubic lattice. 

A lattice polygon has length n if it is composed of n edges and n vertices. The 
number of lattice polygons of length n is the number of distinct polygons of length n, 
denoted by p„. The function p„ is the most basic combinatorial quantity associated 
with lattice polygons, and logp„ is a measure of the entropy of the lattice polygon at 
length n. 

Determining p„ in regular lattices is an old and difficult combinatorial problem 
[T6] . Observe that P2n+i = for n S N in the SC lattice, and it known that the 
growth constant n defined by the limit 

hm pi/" = M > (1) 

n— >-oo 

exists and is finite in the SC lattice [TB] if the hmit is taken through even vahies of 



n. This result can be extended to other lattices, including the FCC and the BCC 
lattices, using the same basic approach in reference [16j (and by taking limits through 
even numbers in the BCC). In the hexagonal lattice it is known that /i = v2 + \/2 

m- 

In three dimensional lattices polygons are models of ring polymers. Knotted 
polygons are similarly a model of knotted ring polymers, see for example reference 
|13j on the importance of topology in the chemistry of ring polymers, and ^9] on the 
occurrence of knotted conformations in DNA. 

1.1. Knotted Polygons 

Let S^ be a circle. An embedding of S^ into M.^ is an injection / : 5^ -> E"^. We say 
that / is tame if it contains no singular points, and a tame embedding is piecewise 
linear and finite if the image of / is the union of line segments of finite length in M.^ . 
A tame piecewise linear embedding of S^ into M.^ is is also called a polygon. 

Tame embeddings S^ into M"^ are tame knots, and the set of polygons compose a 
class of piecewise linear knots denoted by ICp. If the class of all lattice polygons (for 
example, in a lattice L) is denoted by V, then V C JCp so that each lattice polygon 
is also a tame and piecewise linear knot in R^ . This defines the knot type K of every 
polygon in a unique way. In particular, two polygons in V are equivalent as knots if 
and only if they are ambient isotopic as tame knots in ICp. 

Define Pn{K) to be the number of lattice polygons in L, of length n and knot type 
K, counted modulo equivalence under translations in L. Then pn{K) is the number of 
unrooted lattice polygons of length n and knot type K. Observe that Pn{K) = if n is 
odd, and hence, consider Pn{K) to be a function on even numbers; Pn{K) : 2N ^ N. 

It follows that p„(Oi) = if n < 4 and Pi{Qi) = 3 in the SC lattice where Oi is 
the unknot (the simplest knot type). \i K ^ Qi is not the unknot, then in the SC 
lattice it is known that Pn{K) = if n < 24 and that P24{K) > J^. In particular, 
P24{3i) = 3328 [43] while Pn{K) ^ ii K :^ Oi or K ^ 3i. 

It is known that 

limsup [pn{K)]^'" ^JIj^ < ^1 (2) 

in the SC lattice; see reference [1B|. li K = Oi is the unknot, then it is known that 

lim [p„(Oi)]^/" = ^0 < M (3) 

n— >oo 

and it follows in addition that /iq < J^k < M! see for example |221 I23j . There are 
substantial numerical evidence in the literature that /Iq = Jij^ for all knot types K 
(see reference [53] for a review, and references [HI [Ml [35] for more on this). Overall, 
these results are strong evidence that the asymptotic behaviour of p„ [K) is given to 
leading order to 

p^{K) ~ CKn"°+^--VS, (4) 

where Nk is the number of prime components of knot type K, and ao is the 
entropic exponent which is independent of knot type. The amplitude Ck is dependent 
on the knot type K. In particular, simulations show that the amplitude ratio 
[pn[K) / pn{L)] — >■ [Ck/Cl] 7^ if Nk = Nl [S]; this strongly supports the proposed 
scaling in equation Q. 

Growth constants for knotted polygons in the FCC and BCC in equations ^ 
and ([3| have not been examined in the literature, but there is general agreement that 



Figure 2. Concatenating polygons in tiie SC lattice. The top edge of the polygon 
on the left is defined at that edge with lexicographic most midpoint, and the 
bottom edge of the polygon on the right as that edge with lexicographic least 
midpoint. By translating, and rotating the polygon on the right until its bottom 
edge is parallel to the top edge of the polygon on the left, and translated one step 
in the X-direction, the two polygons can be concatenated into a single polygon by 
inserting the dotted polygon of length four between the two, as illustrated, and 
then deleting edges which are doubled up. If the polygon on the left has length 
n and knot type K, then it can be chosen in p„{K) ways, and if the polygon on 
the right has length m and knot type L, then it can be chosen in pm{L)/2 ways, 
since its bottom edge much be parallel to the top edge of the polygon on the left. 
This shows that p„{K)pm,{L) < 2p„+m{Ki^L), since the concatenated polygon 
has length n + m and knot type the connected sum K^L of the the knot types 
K and L. This construction generalises in the obvious way to the FCC and BCC 
lattices. 



the methods of proof in the SC lattice will demonstrate these same relations in the 
FCC and BCC. In particular, by concatenating SC lattice polygons as schematically 
illustrated in figure [2J it follows that 

p,,{K)p^(L) < 2p^+„,{K#L) (5) 

where K^L is the connected sum of the knot types K and L. 

Similar results are known in the FCC lattice: One has that Pn{Oi) = if n < 3, 
and P3(0i) — 8. Similarly, p„(3i) = if ?i < 15, while ^15(81) = 64. Observe that 
in the FCC lattice, Pn{K) is a function on N; Pn{K) : N ^- N. That is, there are 
polygons of odd length. 

The construction in figure [2] generalises to the FCC lattice. In this case, the top 
vertex of the polygon is that vertex with lexicographic most coordinates. The top 
vertex t is incident with two edges, and the top edge is that edge with midpoint with 
lexicographic most coordinates. The top edge of a FCC polygon is parallel to one of 
six possible directions, giving six different classes of polygons. One of these classes is 
the most numerous, containing at least Pn{K)/6 polygons and with top edge parallel 
to (say) direction f3, if the polygons has length n and knot type K. 

Similarly, the bottom vertex and bottom edge of a FCC polygon of length m 
and knot type L can be identified, and there is a direction 7 such that the class of 
FCC polygons with bottom edge is parallel to 7 is the most numerous and is at least 
p„(L)/6. 

By choosing a polygon of knot type K, top vertex t and length n with top edge 
parallel to /3, and a second polygon of length m, bottom vertex 6, with bottom edge 
parallel to 7, these polygons can be concatenated similarly to the construction in 
figure [2] by inserting a polygon of length (say) iV + 2 between them. Accounting for 
the number of choices of the polygons on the left and right, and for the change in the 
number of edges, this shows that 

Pn{K)pra{L) < 36 Pn+m+NiK#L) (6) 

in the FCC, where TV is independent of n and m . The polygon of length A^ + 2 is 
inserted to join the top and bottom edges of the respective polygons, since they may 



not be parallel a priori to the concatenation. Some reflection shows that the choice 
A^ = 2 is sufficient in this case. 

The relation in equation (|6| shows that [p„_Ar/36] and [p„_Ar(0i)/36] are 
superniultiplicative functions in the FCC, and this proves existence of the limits 
lim„^oo [Pri] " = /-« and lim„^.oo [Pn(Oi)] = /io in the FCC [19 . In addition, 
with Jlj^ defined in the FCC as in equation ^ , it also follows from equation ^ that 
A*Oi ^ J^K ^ M- That Ji^ < n would follow from a pattern theorem for polygons in 
the FCC (and it is widely expected that the methods in reference [351 ISZ] will prove 
a pattern theorem for polygons in the FCC). 

In the BCC lattice one may verify that p„(Oi) = if n < 4, and P4{0i) — 12. 
Similarly, p„(3i) = if n < 18, while Pi8(3i) = 1584. Observe that in the BCC lattice 
Pn{K) is a function on even numbers; Pn{K) : 2N — > N, similar to the case in the SC 
lattice. 

Finally, arguments similar to the above show that in the BCC lattice there exists 
an N independent of n and m such that 

Pn{K)p^{L) < 16pn+m+N{K#L). (7) 

Thus, in the BCC one similarly expects that lim„_j.oo [p„] = fi and 

lim„_j.oo [Pn(Oi)] " — /io exists in the BCC, and with JIj^ defined in the BCC as 
in equation (pi), it also follows from equation (16]) that /io < Jlx < /*■ Similarly, a 
pattern theorem will show that /I^ < /i. In the BCC one may choose N — 2. 

Generally, these results are consistent with the hypothesis that /I^ = fj,Q in the 
BCC and FCC lattices, while the asymptotic form for Pn{K) in equation ^ is expected 
to apply in these lattices as well. By computing amplitude ratios [Ck/Cl] in reference 
[221 for a selection of knots, strong numerical evidence for equation (|4| in the BCC 
and FCC were obtained. 

1.2. Minimal Length Knots and the Lattice Edge Index 

Given a knot type K there exists an n^ such that pni^{K) > but pn{K) = for 
all n < nx • The number hk is the minimal length of the knot type K in the lattice 
[TUl E7] . For example, ii K — 3^ (a right-handed trefoil knot) then in the SC lattice 
it is known that P2i{3f) = 1664 while Pni^t) = ^^r all n < 24. Thus n^+ = 24 
is the minimal length of (right-handed) trefoils in the SC lattice i45|. Observe that 
713+ = rig- (= nsi), and this is generally true for all knot types. 

Similar results are not available in the BCC and FCC, although numerical 
simulations have shown that n^+ = 18 in the BCC and n^+ = 15 in the FCC 

[smsaiMj- 

The construction in figure [2] shows that 

nK#K <'2nK and uk^l < riK + ul (8) 

in the SC lattice. More generally, observe that for non-negative integers p, 

ukp <pnK- (9) 

This in particular shows that the minimal lattice edge index defined by 



inf = lim = a^ (10) 

p \_ p j p-¥oo I p \ 

exists, and moreover, ax > 4(6^ — 1), where bx is the bridge number of the knot 
type K; see references [17| HSl [23] for details. Since 6if > 2 if iiT 7^ 0, it follows that 



OiK > 4 for non-trivial knots types in the SC lattice. Observe that ao^ — and that 
it is known that 4: < a^+ < 17 pTll^U] , 

In the BCC and FCC lattices one may consult equations ^ and (l7| to see that 
for non-negative integers, 

nKp#Ki < nKp + riKi + N. (11) 

Thus, ukp -I- A^ is a subadditive function of p, and hence 



inf 

p 



TlKP +N 
P 



lim 



aK (12) 



nRp 

. P 

exists [H]. Moreover, as in the SC lattice, one may present arguments similar to those 
in the proof of theorem 2 in reference [H] to see that a^ > 3(6^ ^ 1) in the FCC and 
ocK > 2(6i<- — 1) in the BCC. Hence, if K is not the unknot, then uk > 3 in the FCC 
and aif > 2 in the BCC. 

We shall also work with the total number of distinct knot types K with uk < "., 
denoted by Qn- It is known that Qn = 1 if n < 24 in the SC lattice, and that Qn = 3 
if 24 < n < 30 |3S], also in the SC lattice. Q„ grows exponentially with n. 

1.3. The Entropy of Minimal Length Knotted Polygons 

If n = nx, then pn{K) > for a given knot type. The entropy of the knot type K 
at minimal length is given by log Pn{K) when n — uk^ More generally, the entropy 
of lattice knots of minimal length and knot type K can be studied by defining the 
density of the knot type K at minimal length by 

Vk=Pu^{K). (13) 

Then one may verify that T^g = 3 in the SC lattice, and 7^0 = 12 in the BCC lattice 
while 7^0 = 8 in the FCC lattice. 

It is also known that P^+ = 1664 in the SC lattice |3S]. Since 3i is a chiral knot 
type, it follows that the total number of minimal length lattice polygons of knot type 
3i is given by "Psi = V^+ + V^- = 3328. 

Generally there does not appear to exist simple relations between Vk and Vk^ . 
However, "P^fm should increase exponentially with m, since Uk^ is bounded linearly 
with m if ii' is a non-trivial knot type [21]. Thus, the entropic index per knot 
component of the knot type K can be defined by 

'logVK" 



lim sup 



= 1K- (14) 



Obviously, since 7'0m — 4 for all values of to, it follows that 70 = 0. Also, ^k > for 
all knot types K. Showing that 7^ > for all non-trivial knot types K is an open 
question. 

The collection of Vk minimal length lattice knots are partitioned in symmetry 
(or equivalence) classes by rotations and reflections (which compose the octahedral 
group, which is the symmetry group of the cubic lattices). Since the group has 24 
elements, each symmetry class may contain at most 24 equivalent polygons. The total 
number of symmetry classes of minimal length lattice knots of type K is denoted by 
Sk- For example, in the SC lattice it is known that S^^ — 1 and this class has 3 
minimal length lattice knots of length 4. It has been shown that ^3^ — 142 in the SC, 
of which 137 classes have 24 members each and 5 have 8 members each [JS]. 

II Sometimes, this notion will be abused when we refer to pn{K) as the (lattice) entropy of polygons 
of length n and knot type K. 





Figure 3. A negative crossing (left) and a positive crossing (rigiit) of the 
intersections in a regular projections of a simple closed curved. The signs are 
assigned using a right hand rule. 



1.4- The Mean Absolute Writhe of Minimal Length Knotted Polygons 

The writhe of a closed curve is a geometric measure of its self-entanglement. It is 
defined as follows: The projection of a closed curve in M.^ onto a geometric plane is 
regular if all multiple points in the projection are double points, and if projected arcs 
intersect transversely at each double point. 

Intersections (referred to as "crossings") in a regular projection are signed by the 
use of a right hand rule: The curve is oriented and the sign is assigned as illustrated 
in figure I3] The writhe of the projected curve is the sum of the signed crossings. The 
writhe of the space curve is the average writhe over all possible regular projections of 
the curve. For a lattice polygon uj this is defined by 

Wr{uj)= -^ / Wr{uJ,u) (15) 



47r j„gs2 

where Wr{ui,u) is the writhe of the projection along the unit vector u (which takes 
values in the unit sphere 5'^ - this is called the writhing number of the projection). 
This follows because almost all projections of w are regular. 

The writhe of a closed curve was introduced by Fuller [T3]. It was shown by 
Lacher and Sumners [55] that the writhe of a lattice curve is given by the average of 
the linking number of w with its push-offs uj + eu, for u ^ S^, and e > small. That 
is, 

Wriuj)^^ [ Lk{uj,uj + eu). (16) 

47r J„gs2 

In the SC lattice, this simplifies to the average linking number of w with four of its 
push-offs into non-antipodal octants: 

I 4 

Wr{uj)^-Y,Lk{uJ,0O + U^) (17) 

4=1 

where, for example, one may take m — (0.5,0.5,0.5), U2 — (0.5,-0.5,0.5), U3 = 
(-0.5,0.5,0.5) and M4 = (-0.5,-0.5,0.5). This shows that 4Wr{K) is an integer. 
The average writhe {Wr{K))n of polygons of knot type K and length n is defined 

by 

{Wr{K))„ = ^— Y. ^rH (18) 

\U}\—71 

where the sum is over all polygons of length n and knot type K. If K is an achiral 
knot, then {Wr{K))n — for each value of n [26] . 



The average absolute writhe (|W^r(iir)|)„ of polygons of knot type K and length 
n is defined by 

{\Wr{K)\)r.= ^^ J2 r.HI (19) 

\u}\—n 

where the sum is over all polygons of length n and knot type K. 

The averaged writhe Wk and the average absolute writhe \W\k of lattice knots of 
both minimal length and knot type K are defined as the average and average absolute 
writhe of polygons of knot type K and minimal length: 

The writhe of polygons in the BCC and FCC lattices can also be determined by 
computing linking numbers between polygons and their push-offs 39J. Normally, the 
writhes in these lattices are related to the average writhing numbers of projections of 
the polygons onto planes normal to a set of given vectors. 

The writhe of a polygon in the FCC lattice is normally an irrational number |15) . 
The prescription for determining the writhe of polygons in the FCC lattice can be 
found in reference P^ and is as follows: Put a = 3arcsec3 — tt and /3 — (7r/2 — q:)/3. 
Then the writhe of a polygon uj is given by 

M^,.(w) = ^ (a^l^r(w,u»)+/3^W^r(^,^'»)) (21) 

where the vectors Ui are defined by Ui — (±3/\/22, ±3/a/22, 2/\/22) for all possible 
choices of the signs, and the vectors Vi are defined by Vi = (±a/5/a/6, ±1/^30, 2/\/30), 
(±1/^30, ±^/5/^/6,2/^/30), (±l/v^, ±1/^38, 6/7^8) again for aU possible choices 
of the signs. The writhing number Wriio^Ui) of to is defined as before as the sum of 
the signed crossings in the projected cj on a plane normal to Ui. 

In the BCC lattice the writhe of a polygon lo can be computed by 

1 ^^ 

^-H = Y^E^'-(^'^') (22) 



where the vectors Ui are defined by u^ = (±1/V10, 3/VlO, 0), (±1/V10, 0, 3/VlO), 
(0, ±1/^10,3/^10), (0,3/N/lO,±l/N/n)), (3/^l0,±l/VT0,0), (3/%/I0,0,±1/n/10), 
for all possible choices of the signs. 

By appealing to the Calugareanu and White formula Lk —T-uj + Wr [3 \5U\ for a 
ribbon, one can compute Wr{uJ, Ui) by creating a ribbon (w, w + eui) with boundaries 
UJ and UJ + eui (this is a push-off of w by e in the (constant) direction oi Ui). Since the 
twist of this ribbon is zero, one has that Wr{uj,Ui) = Wr{uj-\-eui,Ui) = Lk{uj,uj-\-eUi), 
and the writhe can be computed by the linking number of the knot Wr{uJ^ Ui) and its 
push-off Wr{uj, Ui) + eUi. 



Equation (22 1 shows that l2Wr{uj) is an integer in the BCC lattice. Thus, the 



mean writhe of a finite collections of polygons in the BCC lattice is a rational number. 

1.5. Curvature of Lattice Knots 

The total curvature of an SC lattice polygon is equal to tt/2 times the number of right 
angles between two edges. The average total curvature of minimal length polygons of 
knot type K is denoted in units of 27r by /Ca' (that is, the average total curvature is 



2'itK,k)- Obviously A^Oi = 1 in the SC lattice, since every minimal length miknotted 
polygon of length 4 is a unit square of total curvature 2tt. For other knot types the 
total curvature of a polygon is an integer multiple of 7r/2, and the mean curvature is 
thus a rational number times 27r. Hence, for a knot type K, the average curvature of 
minimal length length polygons of knot type K is given by 

{Ck) = 27tJCk (23) 

where JCk is a rational number. 

Similar definitions hold for polygons in the BCC and FCC. In each case the lattice 
curvature of a polygon is the sum of the complements of angles inscribed between 
successive edges. 

In the FCC the curvature of a polygon is a summation over angles of sizes 0, 7r/3 
and 27r/3. Hence 2'kK,k is a rational number similarly to the case in the SC. This 



gives a similar definition to equation ( 23 1 of Kk for minimal length lattice polygons 
of knot type K. Obviously, /Cqi = 1 in the FCC, since each minimal lattice polygon 
of knot type Oi is an elementary equilateral triangle. 

The situation is somewhat more complex in the BCC lattice. The curvature of a 
polygon is the sum over angles of sizes arccos(l/-\/3), tt — arccos(l/A/3) and 0. This 
shows that the average curvature of minimal length polygons of knot type K is of the 
generic form 

{Ck) = Bk arccos(l/V3) + 2it1Ck (24) 

where Bk and ICk are rational numbers. By examining the 12 minimal length 
unknotted polygons of length 4 in the BCC, one can show that Bq-^ = —2 and 
/Coi = 3/2. 

The minimal lattice curvature Ck (as opposed to the average curvature) of SC 
lattice knots were examined in reference [28jp[| For example, it is known that Co^ — 2'k 
while Csj^ — Gtt in the SC lattice |28j . Bounds on the minimal lattice curvature in the 
SC lattice can also be found in terms of the minimal crossing number Ck or the bridge 
number hK of a knot. In particular, Ck > max { (3 + ^/%~+%Ck) 7I'/4 , iirhK^- These 
bounds are in particular good enough to prove that Cg^^ = 97r. A minimal lattice 
curvature index vk is also proven to exist in reference |28| , in particular 

lim = Vk (25) 

n— fcxD Tl 

exists and Ca'" > uvk- It is known that I'Oi = but that 27r < v^-^ < 37r in the SC 
lattice, and one expect that 2'kJCk" > Ck^ > Svrn in the SC lattice. This shows that 
/Cx" increases at least as fast as n in the SC lattice. For more details, see reference 



2. GAS Sampling of knotted polygons 

Knotted polygons can be sampled by implementing the GAS algorithm (301 . 
The algorithm is implemented using a set of local elementary transitions (called 
"atmospheric moves" [25]) to sample along sequences of polygon conformations. The 
algorithm is a generalisation of the Rosenbluth algorithm 44J , and is an approximate 
enumeration algorithm [Ml [25] . 

^ Observe that the minimal lattice curvature of a lattice knot does not necessarily occur at minimal 
length. 
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Figure 4. BFACF elementary moves on polygons in the cubic lattice. These 
(reversible) moves are of two types: Type I decreases or increases the length of 
the polygon by two edges, while Type II is a neutral move which maintains the 
length of the polygon. A move which increases the length of the polygon is a 
positive move, while negative moves decrease the length of the polygon. 



The GAS algorithm can be implemented in the SC lattice on polygons of 
given knot type K using the BFACF elementary moves pQ [H |S] to implement the 
atmospheric moves [3TJ [32] • These elementary moves are illustrated in figure [s] This 
implementation is irreducible on classes of polygon of fixed knot type [35] . 

The BFACF moves in figurel4]are either positive (increase the length of a polygon) , 
neutral (leave the length unchanged) or negative (decrease the length of a polygon). 
These moves define the atmosphere of a polygon. The collection of possible positive 
moves constitutes the positive atmosphere of the polygon. Similarly, the collection 
of neutral moves composes the neutral atmosphere while the set of negative moves is 
the negative atmosphere of the polygon. The the size of an atmosphere of a polygon 
uj is the number of possible successful elementary moves that can be performed to 
change it into a different conformation. We denote the size of the positive atmosphere 
of a polygon a; by a+(w), of the neutral atmosphere by ao(w), and of the negative 
atmosphere by a_(a;). 

The GAS algorithm is implemented on cubic lattice polygons as follows (for more 
detail, see references [3T1 132] )■ Let wq be a lattice polygon of knot type K, then 
sample along a sequence of polygons {ujo, oji, a;2, . . .) by updating coi to Wi+i using an 
atmospheric move. 

Each atmospheric move is chosen uniformly from the collection of possible moves 
in the atmospheres. That is, if Wj has length £j then the probabilities for positive, 
neutral and negative moves are given by 

Pr(+) ex f3e.a+{ujj), Pr(0) ex ao(wj), andPr(— ) ex a_(wj) (26) 

where the parameters /3i were introduced in order to control the transition probabilities 
in the algorithm. It will be set in the simulation for "flat sampling". That, it will be 
chosen approximately equal to the ratio of average sizes of the positive and negative 



atmospheres of polygons of length i: [Sg 






This choice makes the average 



probability of a positive atmospheric move roughly equal to the probability of a 
negative move at each value of i. 

This sampling produces a sequence (ujj) of states and we assign a weight 



W{uJn) 






(ij 



u^i: 



l) 



(27) 



to the state w„. The GAS algorithm is an approximate enumeration algorithm in the 
sense that the ratio of average weights of polygons of lengths n and m tends to the 



10 



ratio of numbers of such polygons. That is, 

W™ P^iKY ^ ' 

The algorithm was coded using hash-coding such that updates of polygons and polygon 
atmospheres were done in 0(1) CPU time. This implementation was very efficient, 
enabling us to perform billions of iterations on knotted polygons in reasonable real 
time on desk top linux workstations. Minimal length polygons of each knot type were 
sieved from the data stream and hashed in a table to avoid duplicate discoveries. 
The lists of minimal length polygons were analysed separately by counting symmetry 
classes, and computing writhes and curvatures. 

Implementation of GAS sampling in the FCC and BCC lattices proceeds similar to 
the implementation in the SC lattice. It is only required to define suitable atmospheric 
moves analogous to the SC lattice moves in figure |4] and to show that these moves 
are irreducible on classes of FCC or BCC lattice polygons of fixed knot types. 

The BCC lattice has girth four, and local positive, neutral and negative 
atmospheric moves similar to the SC lattice moves in figure [4] can be defined in a 
very natural way. These are illustrated in figure [5J Observe that the conformations 
in this figure are not necessarily planar, in particular because minimal length lattice 
polygons in the BCC lattice are not necessarily planar. This collection of elementary 
moves is irreducible on classes of unrooted lattice polygons of fixed knot type K in 
the BCC lattice ^^. 

In the FCC lattice the generalisation of the BFACF elementary moves is a single 
class of positive atmospheric moves and their inverse, illustrated in figure [6J This 
elementary move (and its inverse) is irreducible on classes of unrooted lattice polygons 
of fixed knot type K in the FCC lattice [52]. The implementation of this elementary 
move using the GAS algorithm is described in references [3T| |32J HH] . 



la: ^ 




lb: y y ^-^ 

% •- •- • • • • • 



lib- y ^^ 

Figure 5. Elementary moves on polygons in the BCC lattice. These (reversible) 
moves are of two types: Type I decreases or increases the length of the polygon 
by two edges, while Type II is a neutral move which maintains the length of the 
polygon. 
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la: 




Figure 6. The Elementary move on polygons in the FCC lattice. This is the 
only class of elementary moves in this lattice, there arc no neutral moves. 



3. Numerical Results 

GAS algorithms for knotted polygons in the SC, BCC and FCC lattices were coded 
and run for polygons of lengths n < M where 500 < M < 700, depending on the knot 
type (the larger values of M were used for more complicated compound knots). In 
each simulation, up to 500 GAS sequences each of length 10^ states were realised with 
the purpose of counting and collecting minimal length polygons. In most cases the 
algorithm efficiently found minimal conformations in short real time, but a few knots 
proved problematic, and in particular compound knots. For example, knots types 
(3f )^ and (4i)^ required weeks of CPU time, while (3^ )^(3j")^ proved to be beyond 
the memory capacity of our computers. 

Generally, our simulations produced lists of symmetry classes of minimal length 
knotted polygons in the three lattices. Our data (lists of minimal length knotted 
polygons) are available at the website in reference [51] . 

3.1. Minimal Knots in the Simple Cubic Lattice 

3.1.1. Minimal Length SC Lattice Knots: The minimal lengths hk of prime knot 
types K are displayed in table IT] We limited our simulations to prime knots up to 
eight crossings. In addition, a few knots with more than eight crossings were included 
in the table, including the first two knots in the knot tables to 12 crossings, as well 
as 942 and 947. The minimal lengths of some compound knots (up to eight crossings), 
as well as compound trefoils up to (3^)^ and figure eights up to (4i)^, were also 
examined, and data are displayed in table [2] 

The results in tables IT] and [2] confirms data previously obtained for minimal knots 
in the simple cubic lattice, see for example (21] and in particular reference [45] for 
extensive results on minimal length knotted SC polygons. 

The number of different knot types with minimal length n^ < n can be estimated 
and grows exponentially with n. In fact, if Qn is the number of different knot types 
with riK < n, then Q„i > Qn ii m > n. Obviously, Q„ < Pn, so that 

Q = limsup Qi/" < hm p^" = fi (29) 

by equation ([T]). 

On the other hand, suppose that A^ > prime knot types (different from 
the unknot) can be tied in polygons of length m (that is Qm > A^)- Then by 
concatenating k polygons of different prime knot types as in figure [2| it follows that 

QNm > Y^k=o ( ^ ) =2^. In other words 

Qn > 2"/™. (30) 

if n = Nm where A^ is the number of non-trivial prime knot types that can be tied 
in a polygon of length n. For example, if to = 28, then A^ = 1, and thus Q28 > 2. 
Taking n — >■ 00 implies that A^ — > cx) as well so that liminf„_i.oo Qn > 1- 
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Prime Knot types 


4 


Oi 




24 


3i 




30 


4i 




34 


5i 




36 


52 




40 


6i, 62, 63 




42 


819 




44 


7i, 73, 74, 77, 820 




46 


72, 75, 76, 821 




48 


83, 87, 942 




50 


81, 82, 84, 85, 85, 8s, 89, 810, 811, 8 


13, 814, 816, 947 


52 


812, 815, 817, 818 




54 


9i 




56 


92 




60 


lOi, IO2 




64 


111 




66 


II2 




70 


12i, 122 





Table 1. Minimal Length of Prime Knots in the SC lattice. 



In other words, I < Q < n- 

Thence, one may estimate Qn , and increasing n in Q„ should give increasingly 
better estimates of Q. In addition, if Qn approaches a limit bigger than one, then 
Q > 1 and the number of different knot types that can be tied in a polygon of length 
n increases exponentially with n. 

By examining the data in tables [T] and [2] one observes that Q30 = 4 so that 
Q « 4^/30 ^ 1.0472 . . .. Increasing n to 40 gives Q40 = 16, so that Q sa 1.07177. . .. 
If n = 50, then Q^q > 74, hence Q ~ 1.08989 . . .. These approximate estimates of Q 
increases systematically, suggesting the estimates are lower bounds, and that Q > 1. 

The number of distinct knot types with nn = n is (5„ — Qn — Qn-i, and since 
Qn > Qn-1 and Qn = Q"+°(") , it follows that g„ = if 71 is odd, and Q„ = g"+°(") 
for even values of n. 

There appears to be several cases of regularity amongst the minimal lengths of 
knot types in tables fll The sequence of (A'', 2)-torus knots with N > 3 (these are the 
knots {3i, 5i, 7i, 9i, Hi}) increases in steps of 10 starting in 24. Similarly, the sequence 
of twist knots {4i, 61, 81, lOi, 12i} increments in 10 starting in 30, as do the sequence 
of twist knots {52,72,92,112}, but starting in 36. The sequence {41,62,82,102,122} 
also increments in 10, starting at 30 as well. A discussion of these patterns can be 
found in reference [21] (see figure 3 therein). There are no proofs that these patterns 
will persist indefinitely. 

In table [2] the estimated minimal lengths uk of a few compounded knots are 
given. These data similarly exhibit some level of regularity. For example, the family 
of compounded positive trefoils {(3]^)"'} increases in steps of 16 starting in 24. From 
these data, one may bound the minimal lattice edge index of positive trefoils (defined 
in equation (10)). In particular, ax < n-ifp/p, and ii K = 3f and p = 6, then 
< 17^. This does not improve on the upper bound given in 
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Compound Knot Types 


40 


(3+)^ (3+)(3r) 


46 


(3^)(4i) 


50 


(5+)(3+), (5+)(3r), (5+)(3r) 


52 


(4i)2, (5+)(3+) 


56 


(3+)3, (3+)2(3r) 


72 


(3+)4 


74 


(4+)3 


88 


(3+)^ 


96 


(4i)^ 


104 


(3+)6 


118 


(4i)^ 


140 


(4i)« 



Table 2. Minimal Length of Compound Knots in the SC lattice. 



references |57j and [2D], but if the increment of 16 persists, then ii p = 10 one would 
obtain a^+ < 16| < 17. Prehminary calculations indicated that finding the minimal 

edge number for (3^")^" would be a difficult simulation, and this was not pursued. At 
this point, the argument illustrated in figure 4 in reference pi) proves that a^+ < 17, 
and the data above suggest that a^^ = 24 + 16?! for n < 6. If this pattern persists, then 
aK would be equal to 16, but there is no firm theoretical argument which validates 
this expectation. 

Similar observations apply to the family of compounded figure eight knots {(4i)"}. 
The minimal edge numbers for n < 6 are displayed in table [2] and increments by 22 
such that a4^ = 30 + 22n for n < 6. This suggest that 04^ = 22, but the best upper 
bound from the data in table pi is 23 1. 



3.1.2. Entropy of minimal lattice knots in the SC lattice: Minimal length lattice knots 
were sieved from the data stream, then classified and stored during the simulations, 
which was allowed to continue until all, or almost all, minimal length lattice were 
discovered. In several cases a simulation was repeated in order to check the results. 
We are very sure of our data if Vk ^ 1000, reasonable certain if 1000 < Vk ^ 10000, 
less certain if 10000 < Vk ^ 100000, and we consider the stated value of Vk to be 
only a lower bound if Vk > 100000 in table [S] 

Data on entropy, lattice writhe and lattice curvature, were collected on prime 
knot types up to eight crossings, and also the knot types 9i, 9i, 942, 947, lOi and IO2. 
The SC lattice data are displayed in table [3] As before, the minimal length of a knot 



type K is denoted by uk, and Vk 



- {K) is the total number of minimal length 



SC lattice knots of length uk- For example, there are 3328 minimal length trefoils (of 
both chiralities) of length 713^ = 24. Since 3i is chiral, 7^3+ — 3328/2 = 1664. 

The unknot has minimal length 4, which is a unit square polygon in a symmetry 
class of 3 members which are equivalent under lattice symmetries. The 3328 minimal 
lattice trefoils are similarly partitioned into 142 symmetry classes, of which 137 classes 
has 24 members and 5 classes has 8 members each. These partitionings into symmetry 
classes are denoted by 3^ for the unknot, and 24^'^^8^ for lattice polygons of knot type 
3i (of both chiralities) or 12^'^^4^ for lattice polygons of (say) right-handed knot type 
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Table 3. Data on prime knot types in the SC Lattice. 
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Entropy per unit length of minimal polygons of knot type K is defined by 

8k = '^^. (31) 

UK 

This is a measure of the tightness of the minimal knot. If £k is small, then there 
are few conformations that the minimal knot can explore, and such a knot is tightly 
embedded in the lattice (and its edges are relatively immobile). If £k^ on the other 
hand, is large, then there are a relatively large conformational space which the edges 
may explore, and such a knot type is said to be loosely embedded. 

The unknot has foi = (log3)/4 = 0.27467. . ., which will be small compared to 
other knot types, and is thus tightly embedded. 

The entropy per unit length of the (right-handed) trefoils is £3+ = (log 1664)/24 = 
0.3090..., and it appears that the edges in these tight embeddings are similarly 
constrained to those in the unknot. Edges in the (achiral) knot 4i has £4^ = 
(log3648)/30 = 0.2733. . . < £3+, and are more constrained than those in the trefoil. 
Similarly, for five crossing knots one finds that £=+ = 0.2386 . . . while £r+ — 0.3223 . . .. 

The entropy per unit length seems to converge in families of knot 
types. For example in the (iV, 2)-torus knot family {3^, 5]'", 7^'", 9^} one gets 
{0.3090,0.2386,0.2214,0.2234} to four digits accuracy. Similarly, the family of twist 
knots {3]^,5j,7j,9j} gives {0.3090,0.3044,0.2616,0.2555}, again to four digits. 
Similar patterns are observed for the families {6;,^, 8;,^, 10^} ({0.2008, 0.1876, 0.2059}), 
and {6^,8^,10^} ({0.2427,0.2147,0.2164}). Further extensions of the estimates of 
Vk for more complicated knots would be necessary to test these patterns, but the 
scope of such simulations are beyond our available computing resources. 

Finally, there are some knots with very low entropy per unit length. These include 
7;^ (0.1246), 7| (0.1007), 7^ (0.1257), 83 (0.08065), 8^ (0.06621), 8^^ (0.09129), 8j*6 
(0.07742) and 8^0 (0.1088). These knots are tightly embedded in the SC lattice in 
their minimal conformations, with very little entropy per edge available. 

The distribution of minimal knotted polygons in symmetry classes in table[3]shows 
that most minimal knotted polygons are not symmetric with respect to elements of 
the octahedral group, and thus fall into classes of 24 distinct polygons. Classes with 
fewer elements, (for example 12 or 8), has symmetric embeddings of the embedded 
polygons. Such symmetric embeddings are the exception rather than the rule in table 
I3J For example, amongst the listed prime knot types in that table, only eight types 
admit to a symmetric embedding. 

3.1.3. The Lattice Writhe and Curvature: The average writhe Y^k, the average 
absolute writhe |yV|_R- and the average curvature K-k (in units of 27r) of minimal 
length polygons are displayed in table [3j The results are given as rational numbers, 
since these numbers can be determined exactly from the data. Observe that the writhe 
of simple cubic lattice polygons are known to be rational numbers (33 [5^1 |31j , hence 
the average over finite sets of polygons will also be rational. In addition, the average 
writhe Wk is non-negative in table [3] since the right handed knot was in each case 
used in the simulation. 

In most cases in table [3] it was observed that \Wk\ ~ |W|k, with the exception 
of some achiral knots, which have \yV\K > while Wk = 0. The average absolute 
writhe was zero in only two cases, namely the unknot and the knot 83. Generally, the 
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Table 4. SC Lattice Trefoils of lengths 24, 26 and 26. 



average and absolute average writhe of achiral knots are not equal, but the unknot 
and 83 are exceptions to this rule. 

It is known that achiral knots have zero average writhe [33]; and so Wk — 
if the knot type K is achiral. For example, Wsj = '^jUi ~ 3.4417... (for right 
handed trefoils), and hence 3i is a chiral knot type. This numerical estimate for Ws^ 
is consistent with the results of simulations done elsewhere |26l EH , and it appears 
that Wsj is only weakly dependent on the length of the polygons. For example, in 
table Wl the average and average absolute writhe of polygons with knot type 3^ and 
lengths 24, 26 and 28 are listed. Observe that while Wsj and \yV\K do change with 
increasing n, it is also so that the change is small, that is, it changes from 3.44170 . . . 
for n = 24 to 3.45971 ... to 3.46848 ... as n increments from n — 2A to n = 28. These 
numerical values are close to the estimates of average writhes made elsewhere in the 
literature for polygons of significant increased length, and the average writhe seems 
to cluster about the estimate 3.44 ... in those simulations [IHl H1H5]. 

Generally, the average and average absolute writhe increases with crossing number 
in table (IS]). However, in each class of knot types of crossing number C > 3 there are 
knot types with small average absolute writhe (and thus with small average writhe). 
For example, amongst the class of knot types on eight crossings, there are achiral knots 
with zero absolute writhe (83), as well as chiral knot types with average absolute writhe 
small compared to the average absolute writhe of (say) 81. For example, the average 
absolute writhe of 818 is 10/37. The obvious question following from this observation 
is on the occurrence of such knot types: Since there are chiral knot types with average 
absolute writhe less than 1 for knots on 4, 6, 7 and 8 crossings in table [3] would such 
chiral knot types exist for all knot types on C > 6 crossings? 

The curvature of a cubic lattice polygon is a multiple of tt/2, and hence the average 
curvature will similarly be a rational number times 2tt: That is, {Ck) = 2t:1Ck where 
{Ck) is the average curvature of minimal length polygons of knot type K and ICk is 
the rational number displayed in the last column of table |3] For example, the average 
curvature of minimal length lattice trefoils is 6 [|§|] tt ~ 6.962747r. 

The variability in K-k is less than that observed for the writhe Wk in classes 
of knot types of given crossing number in table [3] Generally, increasing the crossing 
number increases the minimal length of the knot type, with a similar increase in the 
number of right angles in the polygon. This increase is reflected in the increase of K-k 
with increasing uk- 

The ratio ICk/tt-k stabilizes quickly in families of knot types. For example, for 
{N, 2)-torus knots, this ratio decreases with increasing uk as {0.25, 0.141, 0.142, 0.141} 
as K increases along {3i, 5i, 7i,9i}. Similar patterns can be determined for other 
families of knot types. For example, for the twist knots K — {4i, 61, 81, lOi}, the 
ratio is also stable, but a little bit lower: {0.145,0.136,0.131,0.136}. 

Finally, it was observed before equation ( 25 ) that the minimal curvature of a 
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Table 5. Minimal Length of knot types in the FCC lattice. 



lattice knot in the SC lattice, Ck, can be defined and that Cq-^ = 27r, Cs^ ~ 6tt and 
C947 = 97r. The average curvatures in table [3] exceeds these lower bounds in general, 



with equality only for the unknot: For example, K.^^ = 6||^7r and JCg^^ 



However, in each of these knot types there are realizations of polygons with both 
minimal length and minimal curvature. 



3.2. Minimal Knots in the Face Centered Cubic Lattice 

3.2.1. Minimal Length FCC Lattice Knots: The minimal lengths n^^ of prime knot 
types K in the FCC are displayed in table [5] Prime knots types up to eight crossings 
are included, together with a few knots with nine crossings, as well as the knots lOi 
and IO2. In general the pattern of data in table [5] are similar to the results in the SC 
lattice in table [l] Observe that while the knot type 6* can be tied with 40 edges in 
the SC lattice, in the FCC lattice 61 and 62 can be tied with fewer edges than 63. 
Similarly, the knot 7i can be tied with fewer edges than other seven crossing knots in 
the FCC lattice, but not in the SC lattice. There are other similar minor changes in 
the ordering of the knot types in table [5] compared to the SC lattice data in table [l] 
Similar to the argument in the SC lattice, one may define Qn to be the number 
of different knot types with uk < n in the FCC lattice. It follows that Q„ < Pn, so 
that 



Q^limsupQ^" < lim 



l/n 



M 



(32) 



by equation ([T]). 

By counting the number of distinct knot types with uk < n in tables [5] one may 

l/n 

estimate Q by computing Qn : Observe that Q3 = 1 and Q15 = 2, this shows that 
Q « 1.041 .... By increasing n, one finds that Q^c, > 37, and this gives the estimate 
Q w 1.10550 .... This is larger than the estimate of Q in the SC lattice, and may be 



18 



some evidence that the exponential rate of growth of Q„ in the FCC lattice is strictly 
larger than in the SC lattice: That is, Qfcc > Qsc- 

Similar to the case in the SC, the number of distinct knot types with uk — n 
is Qn = Qn — Qn-i, and since Q„ > Qn~i and Qn — Q"'+°("-), it follows that 
Q„ = Q"+o("). 

There are several cases of (semi)-regularity amongst the minimal lengths of knot 
types in tableslsj (TV, 2)-torus knots with 3 < iV < 5 (these are the knots {3i, 5i, 7i}) 
have increases in steps of 7 starting in 15. This pattern, however, fails for the 
next member in this sequence, since ng^ = 35, an increment of 6 from 7i. Similar 
observations are true of the sequence of twist knots. The sequence {4i, 61,81} has 
increments of 7 starting in 20, but this breaks down for lOi, which increments by 6 
over nsi- The first three members of the sequence of twist knots {52, 72, 92} similarly 
have increments in steps of 7, and if the patterns above applies in this case as well, 
then this should break down as well. Observe that these results are different from the 
results in the SC lattice. In that case, the patterns persisted for the knots examined, 
but in the FCC lattice the patterns break down fairly quickly. 

3.2.2. Entropy of minimal lattice knots in the FCC Lattice: Data on entropy on 
minimal length polygons were collected of FCC lattice polygons with prime knot 
types up to eight crossings, and also knot the knots 9i, 9i, 942, 947, lOi and IO2. The 
results are displayed in table |6] The minimal length of a knot type K is denoted by 
UK, and Vk — Cn^iK) is the total number of minimal length FCC lattice knots of 
length UK. For example, there are 64 minimal length trefoils (of both chiralitics) of 
length n3j = 15 in the FCC lattice. Since 3i is chiral, V^+ — 64/2 = 32. 

Each set of minimal length lattice knots are divided into symmetry classes under 
action of the symmetry group of rotations and reflections in the FCC lattice. For 
example, the unknot has minimal length 3 and it is a member of a symmetry class 
of 8 FCC lattice polygons of minimal length which are equivalent under action of the 
symmetry elements of the octahedral group. 

The 64 minimal length FCC lattice trefoils are similarly divided into 4 symmetry 
classes, of which 2 classes have 24 members and 2 classes have 8 members each (which 
are symmetric under action of some of the group elements). This partitioning into 
symmetry classes are denoted by 8^ for the unknot, and 24^8^ for the trefoil (of both 
chiralities) . 

Similar to the case for the SC lattice, the reliability of the data in table[6]decreases 
with increasing values of Vk- We are very certain of our data if Vk ^ 1000, reasonable 
certain if 1000 <Vk < 10000, less certain if 10000 <Vk < 100000, and we consider 
the stated value of Vk to be only a lower bound if Vk ^ 100000 in table [6J 

The entropy per unit length of minimal polygons of knot type K is similarly 



defined in this lattice in equation (31). The unknot has relative large entropy: 
foi = (log8)/3 = 0.693147.... 

The entropy per unit length of the (right-handed) trefoil is £3+ = (log32)/15 = 
0.2310 . . ., which is smaller than the entropy of this knot type in the SC. This implies 
that there are fewer conformations per edge, and the knot may be considered to be 
more tightly embedded. 

The entropy per unit length of the (achiral) knot 4i is £4-^ — (log2796)/20 ~ 
0.3968 . . ., and is less than the trefoil (however, in the SC lattice f (3^) > £(4i)). For 
five crossing knots one finds that £r+ = 0.1760 . . . while £,-+ = 0.2587. . .; these are 
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Table 6. Data on prime knot types in the FCC Lattice. 
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Table 7. Data on trefoils of lengths 15, 16 and 17 in the FCC lattice. 

related similarly to the results in the SC lattice. 

The entropy per unit length in the family of (A^, 2)-torus knots {3^ , 5^, 7^ ,9f } 
changes as {0.2310, 0.1760, 0.2628, 0.1304} to four digits accuracy. These results do not 
show the regularity observed in the SC: While the results for {3^,5^,9]''} decreases 
in sequence, the result for 7^ seems to be unrelated. 

The family of twist knots {3^, 5^, 7^, 9^} gives {0.2310,0.2587,0.2544,0.3149}, 
again to four digits, and this case the knot 9^ seems to have a value higher 
than expected. Similar observations can be made for the families {6]'", 8]^, 10]'"} 
({0.3392,0.3623,0.2642}), and {6^,8^,10^} ({0.2901,0.3226,0.2101}). Further 
extensions of the estimates of Vk for more complicated knots would be necessary 
to determine if any of these sequences approach a limiting value. 

Finally, there are some knots with very low entropy per unit length. These include 
7| (0.1290), 9;;^ (0.1106), and ^t^ (0.1210). These knots are tightly embedded in the 
FCC lattice in their minimal conformations, with very little entropy per edge available. 

The distribution of minimal knotted polygons in symmetry classes in tablelHlshows 
that most minimal knotted polygons are not symmetric with respect to elements of the 
octahedral group, and thus fall into classes of 24 distinct polygons. Classes with fewer 
elements, (for example 12 or 8), contains symmetric embeddings of the embedded 
polygons. Such symmetric embeddings are the exception rather than the rule in table 
|6| This is similar to the observations made in the SC lattice. 



3.2.3. The Lattice Writhe and Curvature in the FCC Lattice: The average writhe 
Wk, the average absolute writhe \yV\K and the average curvature JCk (in units of 27r) 
of minimal length FCC lattice polygons are displayed in table [6) The results for the 
average writhe are given in floating point numbers since these are irrational numbers 



in the FCC lattice, as seen for example from equation (21) 



The lattice curvature of a given FCC lattice polygon, on the other hand, is a 
multiple of 7r/3, and thus 27r/Cx, where ICk is average curvature, is a rational number. 
In table |6] the average curvature K-k is given in units of 27r, so that the exact values of 
this average quantity can be given as a rational number. For example, one infers from 
table |6] that the average curvature of the unknot is 27r, while the average curvature of 
3i is 2f(27r) = 5^77. 

Similar to the results in the SC lattice, the absolute average and average absolute 
writhes in table [3] are equal, except for achiral knots. This pattern may break down 
eventually, but persists for the knots we considered. In the case of achiral knots one 
has, as for the SC lattice, Wk = while \W\k > 0. Observe that the average absolute 
writhe of 83 is positive in the FCC, but it is zero in the SC. 

The average writhe at minimal length of 3^ is 3.3245 ... in the FCC, while it is 
slightly larger in the SC, namely 3.4417. Increasing the value of n from 15 to 16 and 
17 in the FCC lattice and measuring the average writhe gives the results in table [Tj 
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which shows that the average writhe increases slowly with n. However, the average 
writhe remains, as in the SC lattice, quite insensitive to n. 

Generally, the average and average absolute writhe increases with crossing number 
in table[6j However, in each class of knot types of crossing number C > 3 there are knot 
types with small average absolute writhe (and thus with small average writhe). For 
example, amongst the class of knot types on eight crossings, there are achiral knots 
with small absolute writhe (83), as well as chiral knot types with average absolute 
writhe small compared to the average absolute writhe of (say) 81. For example, the 
average absolute writhe of 89, 817 and 818 are small compared to other eight crossing 
knots (except 83). 

While the average writhe is known not to be rational in the FCC, it is 
nevertheless interesting to observe that the average writhe of 816 is almost exactly 
8/3 (it is approximately 41472.00002289/15552 = 8.0000000044/3). Similarly, the 
average writhe of the figure eight knot is very close to 13/200 (it is approximately 
12.99108/200). 

The average curvature ICk tends to increase consistently with uk and with 
crossing number of K. The ratio ICx/nK stabilizes quickly in families of knot 
types. For example, for {N, 2)-torus knots, this ratio decreases with increasing uk 
as {0.183,0.159,0.169,0.164} as K increases along {3i,5i, 7i,9i}. These estimates 
are slightly larger than the similar estimates in the SC lattice. Similar patterns can 
be determined for other families of knot types. For example, for the twist knots 
K = {4i, 61,81, lOi}, the ratio is also stable and close in value to the twist knot 
results: {0.163, 0.162, 0.165, 0.164}. 

Finally, the average curvature of the trefoil in the FCC is 5i7r and this is less than 
the lower bound 6tt of the minimal curvature of a trefoil in the SC lattice |25]. The 
minimal curvature of 947 at minimal length in the SC lattice is 97r [21] , but in the FCC 
lattice our data show no FCC polygons of knot type 947 and minimal length n = 35 
has curvature less than IO^tt. In other words, there is no realisation of a polygon of 
knot type 947 in the FCC at minimal length tik = 35 with minimal curvature 97r. The 
average curvature of minimal length FCC lattice knots of type 947 is still larger than 
the these lower bounds, namely IO^tt. 

3.3. Minimal Knots in the Body Centered Cubic Lattice 

3.3.1. Minimal Length BCC Lattice Knots: The minimal lengths uk of prime knot 
types K in the BCC are displayed in table [8j We again included prime knot types up 
to eight crossings, together with a few knots with nine crossings, as well as the knots 
lOi and IO2. 

In general the pattern of data in table [8] is similar to the results in the SC and 
FCC lattices in tables [l] and [8] The spectrum of knots corresponds well up to five 
crossings, but again at six crossings some differences appear. For example, in the 
BCC lattice one observes that rig^ < ngj and ng^ < uq^, in contrast with the patterns 
observed in the SC and FCC lattices. 

The rate of increase in the number of knot types of minimal length uk < n in the 
BCC lattice may be analysed in the same way as in the SC or FCC lattice. Similar to 
the argument in the SC lattice, one may define Qn to be the number of different knot 
types with Ufc < n in the FCC lattice. It follows that Q„ < Pn, so that 

Q = limsup Qi/" < hm p]/" = /i (33) 
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Table 8. Minimal Length of Knots types in the BCC lattice. 

by equation (IT]). 

By counting the number of distinct knot types with jik < n in table [8] one may 
estimate Q: Observe that (54 = 1 and Qis = 2, this shows that Q sa 2^/^^ — 1.035 . . .. 
By increasing n while counting knot types to estimate Qn, one finds that Q^g > 35, 
and this gives the lower bound Q « 35^^'"' — 1.0929.... This is larger than the 
lower bound on Q in the SC lattice, and may again be taken as evidence that Qn 
is exponentially small in the SC lattice when compared to the BCC lattice. That is 

Qsc < Qbcc- 

Similar to the case in the SC lattice, the number of distinct knot types with 
riK = n is Q^ = Qn — Qn-i, and since (5„ > Q„_i and Q„ — Q"+o(")^ it follows that 
Qn = Q"-+°(") for even values of n (note that Q„ = of n is odd, since the BCC is a 
bipartite lattice). 

There are several cases of (semi)-regularity amongst the minimal lengths of knot 
types in tables l5J (A^, 2)-torus knots (these are the knots {3i, 5i, 7i, 9i}) increase in 
steps of 6 or 8 starting in 18. The increments are {8,6,8} in this particular case, 
and there are no indications that this will be repeating, or whether it will persist 
at all. Similar observations are true of the sequence of twist knots. The sequence 
{4i, 61, 81, lOi} seems to have increments of 8 starting in 20, but this breaks down for 
lOi, which increments by 6 over 81. Similar observations can be made for the sequence 
of twist knots {52, 72, 92}. 

3.3.2. Entropy of minimal lattice knots in the BCC lattice: Data on entropy of 
minimal length polygons in the BCC lattice are displayed in table |9] The minimal 



length of a knot type K is denoted by uk , and Vk 



- (K) is the total number of 



minimal length BCC lattice knots of length uk- For example, there are 1584 minimal 
length trefoils (of both chiralities) of length 71,3^ = 18 in the FCC lattice. Since 3i is 
chiral, P3+ = 1584/2 = 792. 

Each set of minimal length lattice knots are divided into symmetry classes under 
the symmetry group of rotations and reflections in the BCC lattice. For example, the 
unknot has minimal length 4 and there are two symmetry classes, each consisting of 



23 



Knot 






Body Centered Cubic Latt 


ice 






TlK 


Vk 


Sk 


Wk 


Mk 


Bk,ICk 


Oi 


4 


12 


2 


6^ 








9 3 
~2, 2 


3i 


18 


1584 


66 


2466 


40 
■^99 


o40 
■^99 


11 J_ 21 
-'--^33' 22 


4i 


20 


12 


2 


62 








16, 


5i 


26 


14832 


618 


246I8 


6 ^^ 

° 1854 


6 **^ 
° 1864 


10^^ 177 
-'-^ 103 ' 206 


52 


26 


4872 


203 


24203 


4 129 
^203 


4I29 
203 


1 7 23 164 
-'- ' 203 ' 203 


61 


28 


72 


4 


242122 


H 


H 


24, 


62 


30 


8256 


344 


24344 


9 30 
•^43 


2i 


9n21 35 
2U43' 43 


63 


30 


3312 


138 


24138 





4 

69 


in 35 56 

-'-^Bg- 69 


7i 


32 


1464 


61 


2461 


9 


9 


24S,1 


72 


32 


24 


1 


241 


6 


6 


28, 


73 


34 


22488 


937 


24937 


T 919 
'2811 


^7 919 
'2811 


25oq7 745 


74 


34 


11208 


468 


24*66122 


r 464 
^467 


1-464 
^467 


9/394 340 
^^467> 467 


75 


34 


8784 


366 


24366 


7 196 
' 549 


7 196 
' 549 


9924 47 
^^61 ' 61 


76 


32 


48 


2 


242 


H 


H 


26, 


77 


32 


24 


1 


241 


3 


2 
3 


24, 


81 


36 


744 


32 


2430122 


2f 


2| 


32, 


82 


38 


118080 


4920 


244920 


^m 


r 782 
-' 1845 


90 153 431 

^"410' 492 


83 


36 


108 


6 


24*62 








32, 


84 


38 


93984 


3916 


243916 


1 4715 
-'- 11748 


1 4715 
^ 11748 


97 955 3849 
2' 195S' 3916 


85 


38 


7392 


318 


242«8i220 


r331 
■J 924 


t; 331 
=^924 


9q11 195 
^^Tl' 308 


86 


38 


9024 


376 


24376 


42S2 


4552 


90 87 117 
2094' T88 


87 


38 


47856 


1994 


241994 


2035 
^2991 


n 2035 
2 2991 


97 59 784 
^' 997> 997 


88 


38 


34656 


1444 


241444 


.112 
^ 361 


,112 
-^361 


9gl77 280 
^"-'361 ' 361 


89 


38 


5712 


238 


24238 





1 
14 


90 55 185 
^^119' 238 


810 


38 


11088 


462 


24462 


9 313 

^462 


9313 
^ 462 


9c6 125 
^°7i 154 


811 


38 


15888 


662 


24662 


4 49 
^ 1986 


4 49 
^ 1986 


97 198 425 
' 331 ' 662 


812 


36 


12 


2 


62 








24, 


813 


38 


17616 


734 


24734 


-. 241 

^ 734 


-1 241 

^ 734 


9r: 180 561 
^367' 734 


814 


38 


16944 


706 


24706 


4sfe 


4353 


9f;205 263 
20353> 353 


815 


38 


4272 


180 


24176^24 


S^o 


^m 


9452 71 
^^89' 89 


816 


38 


1056 


44 


2444 


9 29 

^44 


929 

244 


97 15 23 

22 ' 44 


817 


38 


912 


38 


2438 





114 


9/1 9 31 
24t9, 38 


818 


40 


8820 


384 


24354^22466 





94 
735 


9/ 116 1 81 
24735- -^245 


819 


32 


1110 


48 


2445^2261 


Q 102 
^185 


102 
*T86 


90 29 64 
^■^ 185 ' 185 


820 


34 


117096 


4879 


244879 


9 372 

^4879 


2 372 
^ 4879 


in 4575 1 614 
-'-^4879 ' 4879 


821 


34 


696 


30 


2428122 


4ff 


443 
^87 


9023 14 
^^29' 29 


9i 


40 


80928 


3372 


243372 


11 14113 
^^ 20232 


1 I 14113 
^^ 20232 


Q9 3287 1 2719 
^^3372 ' -'- 3372 


92 


40 


13824 


576 


24576 


7 J^ 
' 192 


7J- 
' 192 


29211 i77 
^^288 ' 96 


942 


36 


2736 


114 


24114 


,^9_ 
^ 171 


^ 171 


24^, If 


947 


40 


68208 


2842 


242842 


2033 
•^2842 


9 2033 
^2842 


99 1 q 296 
^^49' 3 1421 


lOi 


42 


288 


12 


2412 


3f 


3| 


38, 


102 


44 


9816 


409 


24409 


Q 136 
*'409 


Q 136 


34 409 > I409 



Table 9. Data on knots in the BCC lattice. 
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6 BCC lattice polygons of length 4 which are equivalent under action of the elements 
of the octahedral group. 

The 1584 minimal length BCC lattice trefoils are similarly divided into 66 
symmetry classes, each with 24 members. This partitioning into symmetry classes are 
denoted by 24^^ (66 equivalence classes of minimal length 18 and with 24 members). 
Similarly, the symmetry classes of the unknot are denoted 6^, namely 2 symmetry 
classes of minimal length unknotted polygons, each class with 6 members equivalent 
under reflections and rotations of the octahedral group. 

Similar to the case for the SC and FCC lattices, the reliability of the data in table 
loldecreases with increasing values oiVx- We are very certain of our data HVk ^ 1000, 
reasonable certain if 1000 < Vk < 10000, less certain if 10000 < Vk < 100000, and 
we consider the stated value of Vk to be only a lower bound if Vk ^ 100000. 

The entropy per unit length of minimal polygons of knot type K is similarly 



defined in this lattice as in equation (31). The unknot has relative large entropy 



foi = (log 12)/4 = 0.621226 . . ., compared to the entropy of the minimal length unknot 
in the SC lattice. 

The entropy per unit length of the (right-handed) trefoil is £^+ = (log792)/18 — 
0.37080..., which is smaller than the entropy of this knot type in the SC lattice. 
This implies that there are fewer conformations per unit length, and the knot may be 
considered to be more tightly embedded. 

The entropy per unit length of the (achiral) knot 4i is £4^ = (logl2)/20 = 
0.12424. . ., and is very small compared to the values obtained in the SC and FCC 
lattices. In contrast with the FCC, the relationship between the knot types 3^ 
and 4i in the BCC lattice is similar to the relationship obtained in the SC lattice, 
f (3]^) > f (4i). Five crossing knots in the BCC lattice have relative large entropies. 
One finds that £,+ = 0.34274 . . . while £.+ = 0.29992 .... 

The entropy per unit length in the family of (A^, 2)-torus knots {3^ ,5^ ,7^ ,9^} 
changes as {0.3708,0.3427,0.2061,0.2652} to four digits accuracy. These resuhs do 
not show the regularity observed in the SC lattice: While the results for {S^ ,5^ ,7i} 
decreases in sequence, the result for 9^ seems to buck this trend. 

The family of twist knots {3j*', 5^, 7^, 9^} gives {0.3708,0.3000,0.0777,0.2210}, 
again to four digits, and this case the knot 72 seems to have a value lower 
than expected. Similar observations can be made for the families {6^, 8^", lOJ"} 
({0.1280,0.1644,0.1183}), and {6^,8^,10^} ({0.2775,0.2891,0.1932}). Further 
extensions of the estimates of Vk for more complicated knots would be necessary 
to determine if any of these sequences approach a limiting value. 

Finally, there are some knots with very low entropy per unit length. These include 
4i (0.1242), 6t (0.1280), 7^ (0.0777), 7J (0.0993), 7^ (0.0777), and 8j^2 (0.0498). 
These knots arc tightly embedded in the BCC lattice in their minimal conformations, 
with very little entropy per edge available. 

The distribution of minimal length knotted polygons in symmetry classes in table 
|6] shows that most minimal knotted polygons are not symmetric with respect to 
elements of the octahedral group, and thus fall into classes of 24 distinct polygons. 
Classes with fewer elements, (for example 12 or 8), contains symmetric embeddings of 
the embedded polygons. Such symmetric embeddings are the exception rather than 
the rule in table |9l This is similar to the observations made in the SC and FCC 
lattices. 
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Body Centered Cubic Lattice 




n 


Vk 


Sk 


Wk 


Mk 


^K,^K 


3i 


18 
20 
22 


1583 

236928 

21116472 


66 

9879 
879864 


2466 

2498651214 
24879842^222 


3i 

22457 
■^ 59232 
1050094 
■^ 2639559 


■^99 

, 22457 

■^59232 

q 1050094 

^ 2639559 


1 1 4 21 
1^33' 22 
n8369 1 12133 
y9872> 1 19744 
n 41747 n 116513 
^879853 ' " 879853 



Table 10. Data on trefoils of lengths 18, 20 and 22 in the BCC lattice. 

3.3.3. The Lattice Writhe and Curvature: The average writhe Wk, the average 
absolute writhe \W\k and the average curvature JCk (in units of 27r) of minimal length 
BCC lattice polygons are displayed in table [9] The writhe Wr{u}) of a BCC lattice 
polygon w is a rational number (since 12Wr{(j~>) is an integer) as shown in equation 
(22 1. Thus, the average writhe and average absolute writhe of minimal length BCC 



lattice polygons are listed as rational numbers in table |9] These results are exact 
in those cases where we succeeded in finding all minimal length BCC polygons of a 
particular knot type K. 

The lattice curvature of a given BCC lattice polygon is somewhat more 
complicated. Each BCC lattice polygon w has curvature which may be expressed in the 
form B arccos(l/-\/3) + 27rC, where B and C are rational numbers. Thus, the average 
curvature of minimal BCC lattice polygons of knot type K is given by expressions 
similar to equation ( pi] ), with Bk and K-k rational numbers. In table [9] the values of 
Bk and JCk are given for each knot type, as a pair of rational numbers. For example, 
the average curvature of the unknot is —2 arccos(l/-\/3) + 37r = 7.51414 . . . > 27r. This 
shows that some minimal conformations of the unknot are not planar. 

Similarly, the average curvature of minimal length polygons of knot type 3i is 
given by 11;^ arccos(l/-\/3) + jyTt = 16.6218. . . < 67r. In other words, the average 
curvature of minimal length BCC lattice trefoils is less than Gtt, which is the minimal 
lattice curvature of SC lattice trefoils. In fact, one may check that this average 
curvature is less than 5^7r, which is the average curvature for minimal length FCC 
lattice polygons. In other words, the embedding of lattice trefoils of minimal length 
in the BCC has lower average curvature than either the average curvature in the SC 
or FCC lattices. 

Similar to the results in the SC and FCC lattices, the average and average absolute 
writhes in table |3] are equal, except in the case of achiral knots. If K is an achiral 
knot type, then generally Wk = while \W\k > 0, similar to the results in the SC 
and FCC lattices. Observe that the average absolute writhe of 83 is zero in the BCC 
lattice, as it was in the SC lattice (but it is positive in the FCC lattice). 

The average writhe at minimal length of 3^ is 3||7r in the BCC lattice, which 
is slightly smaller than the result in the SC lattice (3 3^). However, it is still larger 
than the result in the FCC lattice. Increasing the value of n from 18 to 20 and 22 
in the BCC lattice and measuring the average writhe of 3^ gives the results in table 
[TOl which shows that the average writhe decreases slowly with n, in contrast with the 
trend observed in the FCC lattice. However, the average writhe remains, as in the SC 
lattice, quite insensitive to n. 

Generally, the average and average absolute writhe increases with crossing number 
in table (IS]). However, in each class of knot types of crossing number C > 3 there are 
knot types with small average absolute writhe (and thus with small average writhe). 
For example, amongst the class of knot types on eight crossings, there are achiral 
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knots with zero absolute writhe (83 and 812), as well as chiral knot types with average 
absolute writhe small compared to the average absolute writhe of (say) 82. The knot 
types 84, 88, 8g, 813, 817 and 818, amongst knot types on eight crossings, also have 
average absolute writhe less than 2, which is small when compared to other eight 
crossing knots such as 82. 

The average curvature of minimal length BCC lattice knots are given in 
terms of the rational numbers Bk and /C^, as explained above. Both Bk 
and K-K tends to increase with rin in table ^ The ratios [K^/n/f ,/C/f/nx] 
however, may decrease with increasing tik within families of knot types. For 
example, for {N, 2)-torus knots, these ratios decrease with increasing riK as 
{[0.618,0.053], [0.745,0.033], [0.770,0.031], [0.824,0.045]} to three digits accuracy 
along the sequence {3i, 5i, 7i,9i}. Similar patterns can be determined for other 
families of knot types. 

Finally, the average curvature of the trefoil in the BCC lattice is 
11t|j arccos(l/\/3) + ffTT = 16.6218 . . . and this is less than the lower bound 67r of 
the minimal curvature of a trefoil in the SC lattice [3S]. The minimal curvature of 
947 at minimal length in the SC lattice is Qtt, but in the BCC our data show that the 
minimal curvature is 28.255468 . . . < Qtt. In other words, there are minimal length 
conformations of the knot 947 in the BCC lattice with total curvature less than the 
minimal curvature 97r of this knot in the SC lattice pSj . 

4. Conclusions 

Data for compounded lattice knots were significantly harder to collect than for the 
prime knot types. Thus, we collected data in only the SC lattice, and we considered 
our data less secure if compared to the data on prime knot types listed in tables |3) [6] 
and[9j The data for compounded SC lattice knots are presented in table [TT] Included 
are the first few members of sequences ((3^)^) and (4]^) and mixed compound knots 
up to eight crossings, with (3^ )^#(3]^) included. We made an attempt to find all 
minimal knots of type (3^)^#(3]^)^ but ran out of computer resources when 7000000 
symmetry classes were detected. 

Compound knot types in the SC lattice tended to have far larger numbers of 
symmetry classes at minimal length, compared to prime knot types with similar 
minimal length or crossing numbers. We ran our simulations for up to weeks in 
some cases, in an attempt to determine good bounds on the numbers of minimal 
length polygons. As in the case of prime knots, certainty about our data decreases 
with increasing numbers of symmetry classes, from very certain when Vk ^ 1000, to 
reasonably certain when the number exceeds 1000 < Vk S 10000, less certain when 
10000 < Vk ^ 100000, and the stated value o^Vk should be considered a lower bound 
if Vk > 100000. 

That is, the data in table [lljfor the knot types ((3+)"), ((3+)^) and (4f) may 
not be exact for Vk, symmetry classes and estimates of the writhe and curvature. At 
best, those results are lower bounds on the counts, within a few percent of the true 
results. 



The data in table 11 allows us to make rough estimates of 73^ (see equation (14)). 
By taking logarithms of Vs^ , one gets for 73^ the following estimates with increasing 
N: {8.1101, 5.1639, 4.1912, 3.8838, 3.4424}. These values have not settled down and it 
is apparent that simulations with more complex compounded knots will be needed to 
estimate jk- 
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Table 11. Data on knots in the SC Lattice. 



In addition, we can make estimates analogous to 7^ by considering the writhe or 
curvature instead: Define 
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K^ 



N 



and 



Pk — lim sup 
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/c 



XN 



N 



(34) 



then one may attempt to estimate these numbers for the trefoil and figure eight 
knots. (^K can be interpreted as the average absolute writhe per knot component, 
and similarly, /3/f is the average curvature at minimal length per knot component. 

The data for the trefoil give the sequence {3.441, 3.459, 3.464, 3.465, 3.468}. These 
results show that C31 ~ 3.47. The similar analysis for 4i gives {0.217, 0.124, 0.099} and 
it appears that there is a more pronounced dependence on the number of components 
in this case. It is difficult to estimate ^4^ from these results, and we have not ruled 
out the possibility that it may approach zero as the number of components increases 
without bound. 

Repeating the above for /Jg^ gives the estimates {3.481, 2.718, 2.428, 2.360, 2.336} 
so that one may estimate /?3j ss 2.3 in units of 27r. Observe that the bounds on v^-^ 
following equation (25) suggest that (3k > 1 for any knot type K y^ Oi. The estimates 



for 4i are {4.007,3.336,2.960}, so that one cannot yet determine an estimate for /S^-^. 
Overall we have examined the entropic and average geometric properties of 
minimal length lattice knots in the SC, the FCC and the BCC lattices. Our data 
were collected using Monte Carlo algorithms with BFACF-style elementary moves. 
The statistical and average properties of sets of minimal length knotted polygons 
were determined and discussed, and comparisons were made between the results in 
the three lattices. Our results show in particular that the properties of minimal length 
lattice knots are not universal in the three lattices. The spectrum of minimal length 
knot types, the entropy, and the average lattice curvature and lattice writhe shows 
variation in several aspects. For example, the spectra of minimal length knots in tables 
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Table 12. Ranking of minimal length lattice knots of knot types to 
by entropy per unit length in the SC, FCC and FCC lattices. 



crossings 



111 [5] and Is] do not maintain a strict order, but shuffle some knots types up or down 
the table in the different lattices. 

Similar observations can be made with respect to the entropy of minimal length 
knots. For example the entropy of the knot types 5i and 62 are inverted in the BCC 
lattice, compared to the relation they have in the SC and FCC lattices. In table [T2| we 
rank knot types by the entropy per unit length at minimal length. That is, we rank 
the knot types by computing f^ (see equation (31)) ~ the larger the result, the lower 
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the ranking in the table (that is, the higher the knot type is Hsted in the table). The 
rankings in table [12] are shuffled around in each of the three lattices. For example, the 
trefoil knot is ranked at position 1 in the SC lattice, at position 26 in the FCC lattice, 
and at position 2 in the BCC lattice. Other knot types are similarly shuffled. 

In the case of writhe there are also subtle, but interesting differences between the 
three lattices. For example, the average absolute writhe of the knot type 83 is zero 
in the SC and BCC lattices, yet it is not zero in the FCC lattice. Equally interesting 
about the results in the FCC lattice is the fact that the average absolute writhes of 
the knot types 4i and Sig are very nearly very simple fractions (far simpler than the 
number of symmetry classes in each case would suggest), in addition to the fact that 
the average absolute writhe of the knot 4i is identically zero in the BCC lattice (but 
not in the SC and FCC lattices). 

Finally, an analysis of the number of knot types of minimal length n^ < n, 
denoted Qn, suggest that Qn ~ Q"- Our data suggest that Qsc < Qbcc < Qfcc, 
so that that the number of knot types which can be tied in a polygon of n edges 
increases fastest (at an exponential rate) in the FCC lattice, followed by the BCC 
lattice and then the SC lattice. 
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